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Abstract 

The COBE satellite has provided the only comprehensive multi-frequency full-sky 
observations of the microwave sky available today. Assessment of the observations requires 
a detailed likelihood analysis to extract the maximum amount of information present in 
the noisy data. I present a specific method for estimating the CMB anisotropy power 
spectrum independent of any assumptions about the underlying cosmology, and then use 
standard image processing techniques to generate the most revealing corresponding maps 
of the signal. The consistency of the data at the three available frequencies provides 
strong support to the assertion that we are being provided with our first glimpse of the 
last scattering surface. 


1 Introduction 

Operations of the Differential Microwave Radiometer instruments (DMR), the last active ex¬ 
periment on board NASA’s COsmic Background Explorer {COBE) satellite, were terminated in 
December 1993, concluding four years of measurements of anisotropy of the cosmic microwave 
background (CMB) radiation. The final product of the DMR-team work — full sky maps of 
the microwave sky at 31.5, 53, and 90 GHz — were released to the astronomical community in 
January 1996. A generic description of the final DMR data set and brief summary of the results 
of the DMR-team analysis were given in [^. The issues of prime concern in the DMR-team 
work included the following: 
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Figure 1: Power spectra of the DMR (galactic frame, extended Galaxy cut, 3881 pixels) differ¬ 
ence maps, (A-B)/2, at each frequency of observations, for each yearly sky maps (top row, stars 
— year one, squares - year 2, triangles — year 3, diamonds — year 4), for two 2-year sky maps 
(middle row, stars — years 1 and 2, diamonds — years 3 and 4), and for the hnal four-year 
product (bottom row, diamonds). Bottom row plots include the 68, 95, and 99% conhdence 
regions (heavy to light grey) from Monte Carlo simulations of instrumental noise in the sky 
maps. 


1) modeling and removal of identified systematic artifacts from the data (see [|T^]) to produce 
sky maps suitable for cosmological studies, 

2) studies of potential signal contamination by the Galaxy [|TI|, and/or the extragalactic sources 

i, 

3) testing of the hypothesis of primordial origin of the measured CMB anisotropy 0] by evalu¬ 
ation of the frequency dependence of the {6T/T)rms in the sky maps, 

4) testing of the hypothesis of gaussianity of statistics of the measured anisotropies 
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5) evaluation and analysis of the auto- and cross-correlation functions of the sky maps Q, 

6) determination of the angular power spectrum of the CMB anisotropy i, 0, 0]- 
Other analyses of the 4 year data presented so far include P, [0, and Bond and Jaffe in these 
Proceedings. 

In this contribution I will focus on two aspects of the interpretation and visualisation of the 
COBE-DMK data: 1) derivation of the angular power spectrum of the CMB anisotropy using 
a data reduction method which is independent of the cosmological model, and 2) linear filtering 
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Figure 2: Power spectra of the DMR (galactic frame, extended Galaxy cut, 3881 pixels) inverse- 
noise-variance weighted sum maps, at each frequency of observations, for all yearly, 2-year, and 
4-year sky maps (same coding as in Figure 1.) Variance weighting of the A and B sides makes 
the yearly and bi-yearly effective noise levels in the 31GHz map more consistent with one 
another than in the corresponding difference maps (with hfty-hfty noise contributions from the 
A and B sides), and the effective noise level in the sum maps somewhat lower than in the 
corresponding difference maps at all three frequencies. 


of the sky maps. Both these goals are pursued within the mathematical framework for analysis 
of the DMR data which was originally described in [Q, and |^. I will briefly describe the DMR 
data, address the necessary technical points, and proceed to presentation of the results and 
discussion. 


2 Data 

The main result of the COBE-DMR. mission is a set of full sky maps, comprised of two maps (A 
and B sides) at each frequency of 31.5, 53, and 90 GHz, pixelized in both galactic and ecliptic 
coordinates. The angular resolution of the instrument was ~ 7°FWHM, and the sky maps are 
divided into 6144 nearly equal area pixels. 

A record of the sensitivity of each detector and the number of observations per pixel in each 
map are provided with the sky maps. This allows for stochastic modelling of the instrument 
noise, which is essential for a proper interpretation of the results of the data analysis and 











subsequent inference of a parametric description for the cosmological signal. [Indeed, except 
for the dipole anisotropy and the galactic plane emission, the CMB anisotropy signals detected 
in the DMR maps are weak compared to typical noise contamination even at 53 GHz — rms 
signal ~ 35 pK vs. rms noise ~ 80 pK in a single pixel.] 

Since the monitoring of the multichannel instrument performance was very long, DMR is 
to date unparalelled among the CMB anisotropy experiments in terms of noise control. 

Using the customized Fourier tools introduced in 0 for the analysis of the incomplete sky of 
the Galaxy-excised DMR sky maps one can extract from the data the cut-sky power spectra of 
the sum and difference maps at each frequency, which are shown in Figures 1 and 2. Quantities 
plotted there are the rms values of the cut-sky harmonic coefficients (the usual {2i -f 1) degrees 
of freedom per Gmode). 

One should appreciate from Figure 1 the steady decrease with time ( 1/Vt) of the rms 

noise amplitude, the excellent consistency of the yearly performances of the 53 and 90 GHz 
instruments, and very good consistency of the measured noise properties with the Gaussian 
model. 

Figure 2 shows the power spectra of signal and noise. Even the totality of the 31.5 GHz 
data remains substantially noisier than the higher frequency channels, which together with un¬ 
certainties in the galactic foreground contribution renders its interpretation more problematic. 
Visual inspection of the 53 and 90 GHz power spectra provides a sense of strong consistency 
between these independent data sets. 


3 Maximum likelihood estimation of the angular power spectrum of 
CMB anisotropy 

A method of cut-sky Fourier decomposition and maximum likelihood analysis of the DMR 
data to infer the parameters of a theoretical spectrum of GMB anisotropy was presented in [^ . 
The application of this formalism to derive the maximum likelihood estimates of the overall 
amplitude and shape parameters for power law spectrum hts to the 2-yr and 4-yr DMR data 
was given in |0|, and §], respectively. Here I demonstrate how this formalism allows one to 
derive a maximum likelihood estimate of the individual multipoles where 

ST/TiVL) = 

A proper analysis of a relatively complex data set such as DMR and the correct inference of a 
parametrisation of the sky distribution of GMB anisotropy requires the simultaneous treatment 
of several interrelated problems, which include: 

1) incomplete sky coverage and resulting multipole coupling after excision of the galactic plane 
region from the analysed maps; 

2) removal of the cosmologically irrelevant low order, i = 0, and 1 modes, preferably without 
affecting the higher order modes, which do carry cosmological information; 

3) modeling of the nonuniform noise distribution; 

4) identification and removal of possible diffuse foreground emission in the areas outside the 
galactic cut; 

5) execution of an unbiased, or least-biased, parameter inference, preferably using the proba¬ 
bility density which describes the tested model of anisotropy in an exact manner. 

A program of analysis of the DMR data which addressed all of these points was conducted in 
0.0. and 1^. An important feature of this method of analysis is that it employs geometrically 
constructed orthogonalized spherical harmonics, which are independent of any cosmological 
theories that one might endeavour to test with the DMR data. This is radically different from 
the possible Karhunen-Loeve type methods, where one builds the basis of modes for sky map 
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Figure 3: Power spectrum of the CMB anisotropy derived by the maximum likelihood method 
with each coefficient treated as a separate parameter (hence Ai/i = 0). Noise and incomplete 
sky coverage induce correlations, impossible to show on the plot, between the derived conhdence 
intervals on each Ci. One should note that a correct estimator of excess variance in signal over 
noise in the data was applied, hence no negative Ci values were obtained. A low-£ tail of the 
CDM model CMB anisotropy spectrum (independently htted to the same data) is shown as a 
slightly tilted solid line. 


decomposition using a specihc theoretical spectrum of CMB anisotropy (e.g. Bond and Jaffe 
in this volume). 

Let us now establish how the method introduced in P can be used to compute an estimator 
of the full sky angular power spectrum, i.e. the coefficients Ci which can be determined from 
the data. Since we work with the Fourier amplitudes c which are linear combinations of 
pixel temperatures (as they are integrals on the cut-sky of the map products with the cut- 
sky orthogonalized spherical harmonic basis functions), their joint probability distribution is 
Gaussian (see [|| for formal details) 

P(c) dc = - = exp f— ■ (S -I- N)“^ ■ c /2) , (1) 

( 27 r)^/ydet(S + N) V ^ ^ V 

where C = S -|- N, and S and N are covariance matrices of the signal and noise, respectively. 
For the maximum likelihood calculations I use the likelihood function dehned as 


F(A|c) = det[S(A) + N] + ■ [S(A) + N]’^ ■ c. 


( 2 ) 










where the full sky signal variance matrix has the form 


A — diag[a 2 , ■ ■ ■, 0,2^ ^31 ■ ■ ■ ■> ^31 

' -V-' '-V-" 
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+1) times 


(3) 


The maximum likelihood problem as formulated in eqs. (1-3) — seeking the coefficients C 2 - 
given the data vector c — involves the assumption of a Gaussian distribution of both noise 
and CMB anisotropies, and the statistical isotropy of the held of temperature perturbations 
(hence only the i dependence of (la^mP)), but does not imply any relation between Ci-s at 
different values of i — hence we are solving for the power spectrum of CMB anisotropy fully 
indpendently of any specihc cosmological model. 

A maximum likelihood search for the power spectrum that best matches the distribution 
of signal on the sky simply entails hnding a minimum of the function F{A) in the imax — 1- 
dimensional space. In the present case it is greatly facilitated by our ability to write down 
easily the expression for the gradient dF/dai, and implement a standard numerical routine 
for an extremum search with gradient zeroing. The solution to this problem for the DMR 4yr 
galactic frame data is shown in Figure 3. Our power spectrum estimator is really an estimator 
of excess variance in signal over noise in the data. One should note in Fig. 3 the asymmetry of 
the conhdence intervals around the maximum likelihood solutions for the C^-s, and the fact that 
no derived values go negative. If there is no signal content in the data that can be attributed 
to a given Ci, an upper limit is established, as e.g. for i = 16. The signal to noise ratio in 
Fourier space drops to ~ 1 at £ ~ 15 in the DMR 4yr data. It is important to see that nearly 
all multipoles within that range are determined with great signihcance. As for the conhdence 
intervals, the plot can not convey the notion of correlations between the errors on different 
multipoles. These correlations are quite small, and more signihcant for multipoles separated 
by Ai = 2 than otherwise. Because of the nontrivial (non-Gaussian) structure of conhdence 
interval correlations it is difficult to implement an exact method of htting of the theoretical 
model of GMB anisotropy to the DMR data by using the power spectrum, compressed rendition 
of the data. For such applications the accurate linear methods, as e.g. |0, should be preferred. 


4 Linear filtering of the DMR sky maps 

When the statistical properties of the signal, polluted by noise in the data, are known, one can 
apply powerful tools of image processing to improve, oftentimes visually quite dramatically, 
the presentation of the measurements. The best known device for linear hltering of noisy 
data is the celebrated Wiener hlter. Alternatives, usually referred to as sub-optimal hlters, 
do exist, and on occasion provide a preferred means of treating the noisy data. In the case 
of the DMR raw data the signal to noise ratio is pretty poor even in the 4yr data set. The 
Wiener hlter, which practically “annihilates” those components of the data which are assesed 
to be noise dominated, when applied to data of low quality usually renders the impression of 
“over-smoothing”. I choose to apply another device, called the power-equalization hlter, which 
is constructed as follows. 

First let us dehne the measurement, signal, and noise vectors, and their covariance matrices: 

m = s -|- n, (s ■ s^) = S, (n • n"^) = N, (m ■ m^) = S -|- N. ( 4 ) 


Next, dehne the hltered data vector, and the hlter matrix: 

f = F- m = F-(s-|-n). 


(5) 





If the filter is constructed from the requirement that the mean square deviation between the 
hltered and the true signal, e = ((f — s)^ ■ (f — s)) = tr((f — s) ■ (f — s)^), is minimized, one 
obtains the Wiener filter 

Fw = S-(S + N)’i. (6) 

Another possibility for sub-optimal filtering, which renders the power equalization filter, 
arises when F is chosen such that 


(f ■ F) = S, i.e. FpE ■ (S + N) ■ FIe = S. 


(7) 


Equation (7) means that on average over many applications the power-equalization hlter renders 
hltered data, whose statistical distribution matches that of the underlying signal. 

Using Choleski decomposition of the relevant matrices, and the following auxiliary matrices, 
we can construct the power equalization hlter as follows: 


S~' =Ls-i -Lf-i, 


S-i — 


s = r^-i-rs-i. 


(S-hN) ^ - L(s+n)-i • Lfs+N)-i) 
and as a result: 


r'(S+N)-i ~L(S+N)-1’ S-|-n 



L 


T 

(S+N)-i 


rfs+N)-i ■ r(s+N)-i, 


( 8 ) 

(9) 

( 10 ) 


This is an upper triangular matrix, hence hltering of a given mode in the cut-sky data vector c 
does not mix in the information from lower order modes, preserving the information ordering 
imposed by the construction of orthonormal modes (see 0). 

For the presentation purposes of this contribution I constructed separate hlters for the 
DMR4 53 and 90 GHz maps. The tremendous advantage of DMR is clearly visible here. We 
can internally compare two sets of decent quality observations at 53 and 90 GHz, and we are 
able to quantify and visualize the consistency with which the instrument reveals the image of 
the last scattering surface. 


5 The first image of last scattering surface 

It should be realized that the GOHF-DMR observations did allow us to make a hrst reliable 
picture of the last scattering surface. Recognition of this fact is somewhat obscured by the 
memory of the picture of the hrst DMR all sky map that was circulated in 1992 - a picture 
which was substantially distorted by noise features. 

Four years after that presentation, the quality of DMR data is now substantially better, 
but it is still somewhat difficult to argue about its validity as an image of the last scattering 
surface, as there is simply no other data at microwave frequencies at the same angular scales to 
allow a direct image to image comparison with DMR. However, since DMR observed the sky 
at three frequencies, and two of these have rendered sky maps of decent quality, an internal 
image comparison is possible. 

What I present here is a result of a fairly long chain of events: separate instruments (with dif¬ 
ferent noise properties) were used to observe the sky (looking at different directions at the time 
of observations!), the data were separately analysed and corrected for the systematic effects, 
the maps were made separately, all the subsequent data reduction was performed separately, 
and in the end different filters were constructed and applied. 

The end product of these operations is visible in two panels on the left side of the color plate 
as the DMR 53 and 90 GHz, power-equalisation filtered sky maps. To facilitate the comparison 


and assessment of consistency of these maps I add and subtract them. This is shown in the 
remaining panels of the color plate in both Mollweide and spherical polar projections. 

It is immediately clear that the 53 and 90 GHz maps contain many common features on 
large angular scales (and, as shown in i, 0, i, they can not be attributed to foreground 
emission from either the Galaxy, or nearby extragalactic objects). The epistemological value of 
this result of the COBE^J^MK mission should not be underestimated, as this is the hrst time 
ever that we are afforded a reliable picture of the most remote regions of the universe. (I shall 
stress that I do not attempt to make an argument about the wonders of image processing, but 
about our truly good luck with the DMR’s splendid performance at two different frequencies!) 

The rms anisotropy in the (53 + 90)/2 PE hltered map is 35/rK, while the rms in the 
(53 — 90)/2 map is about 15/rK. These numbers dehne the accuracy to which COBE-DMK has 
revealed to us a most elusive astronomical realm, and provided a hrst direct glimpse at the 
universe in its embryonic state. 
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